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Abstract 

It is demonstrated experimentally that strongly-coupled plasma exhibits solid superheating. A 
2D suspension of microspheres in dusty plasma, initially self-organized in a solid lattice, was heated 
and then cooled rapidly by turning laser heating on and off. Particles were tracked using video 
microscopy, allowing atomistic-scale observation during melting and solidification. During rapid 
heating, the suspension remained in a solid structure at temperatures above the melting point, 
demonstrating solid superheating. Hysteresis diagrams did not indicate liquid supercooling in this 
2D system. 

PACS numbers: 52.27.Lw, 52.27.Gr, 68.35.Rh, 64.70.D- 
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Strongly-coupled plasma is a collection of free charged particles where the Coulomb in- 
teraction with nearest neighbors is so strong that particles do not easily move past one 
another jl|. Plasma can become strongly coupled due to high density as in neutron stars 2|, 

4|. 



n 



low temperature as in pure ion plasma [3], or high particle charge as in dusty plasma 



Dusty plasma is partially ionized gas containing micron-size particles of solid matter. Dusty 



plasmas have been used in the study of phase transitions 



port 
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13| . and liquid microstructure 
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9(], waves [lOj, trans- 



Materials like water can exist as superheated solid [l5| or supercoo 



the melting point |l^. Observing solid superheating was once t 
but it is now practical due to new instrumentation for heating 
samples 



ed liquid [l6|. These 



are, respectively, a solid at temperatures above the melting point [17|], and a liquid below 

^nought to be impossible [l9| , 
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15|, 



20[ | or fabricating special 



21. 



We find that the literature for solid superheating lacks experiments with atomistic-scale 
observation. Here, the term "atomistic-scale" indicates that molecules or their equivalent 
are imaged or tracked individually. Most solid superheating experiments use external mea- 



surements like diffraction in metals 
measurements for Abrikosov vortices 
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22l | or optical absorption in ice [15[ , or electrical 
23|. In contrast to the experimental literature, theory 



for solid superheating includes simulations that track individual molecules [l7| . Experiments 



with colloidal suspensions include direct imaging of particles in supercooled liquids 
but apparently not superheated solids. 
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25|, 



The literature for solid superheating also lacks experiments with strongly-coiipled plasma. 
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Experiments with strongly-coupled plasma have demonstrated solid and liquid p 
behavior, and recently supercooled liquid as well but not superheated sohd. 

Liquid supercooling, unlike solid superheating, is easily achieved in many three- 
dimensional (3D) systems, but it is an open question whether supercooling ever occurs 



in one-component 2D systems 251]. Experiments are needed to answer this question. Candi- 
date s yste ms for 2D experiments include electrons on a liquid helium surface 26|], granular 
fluids |27|, colloids [251], and dusty plasmas 

Here, we seek answers to three questions. First, can strongly-coupled plasmas exhibit 
solid superheating? Second, can solid superheating experiments be performed using direct 
imaging of particles? Third, does our one-component 2D system exhibit liquid supercooling? 

We report experiments with a 2D suspension of particles in a dusty plasma, which is a 
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kind of strongly-coupled plasma. Highly-charged particles, which are polymer microspheres, 
are immersed in partially ionized argon gas. Electrons and positive ions are collected by 
a particle, giving it a large negative electric charge. In a plane perpendicular to ion flow, 
particles interact through a repulsive Yukawa potential U{r) = (5^(47reor)~^exp(— r/A/j) [28i] . 

Our particles experience multiple forces, the largest arising from gravity, electric fields, 
gas friction, and laser radiation pressure. The apparatus j4| provides a plasma with a sheath 
above a lower horizontal electrode. This sheath has electric fields that levitate and confine 
charged particles, so that they are suspended as a single layer. Particles have a diameter 
4.83 ± 0.08 fim [29] and mass m = 8.93 X 10"^^ kg. To partially ionize 7 mTorr argon 
gas, we used radio-frequency power at 13.56 MHz, with an amplitude of 97 V peak-to-peak. 
Particles experience gas drag with a coefficient of 2.1 s~^ [29] when they move. 

As in colloidal suspensions, our particles can self-organize in a crystal. Unlike colloids, 
however, our particles are underdamped, and they can be heated without heating the gas 
or ions. In our experiment, particle motion was essentially 2D, with negligible out-of-plane 
displacements and no buckling of the particle layer. 

Video microscopy allows imaging this 2D suspension at an atomistic scale, so that we 
can track particles and measure their individual positions and velocities in each video frame. 
Viewing from above, we recorded a movie [sO^ of 5575 frames at 55 frames per second with a 
total field-of-view (FOV) of 34.2 x 25.6 mm^. We analyzed data in a 30.7 x 22.2 mm^ sample 
region in the center of FOV, which included about 1000 of the > 5000 particles in the sus- 
pension. The particle spacing was characterized by a Wigner-Seitz radius {4] of 0.45 mm. In 
each frame, we measured positions of particles and tracked their motion. The particle posi- 
tions were used for three structural indicators, described below. For each video frame, the 2D 
particle velocities Vi were used to calculate the temperature T = {J2f=i ^{vi — vy /2)/NkB, 
where N is the number of particles analyzed, and v is the center-of-mass velocity. This 
kinetic temperature is different from the temperatures of the other constituents including 
the neutral gas, electrons, ions and the polymer material of the particles themselves. Our 
velocity distribution function contained some non-Maxwellian features, as in 311, including 
a peak at = 5 (mm/s)^ as in Fig. 4(c) of [31]. Particle velocities were also used in the 
wave-spectra analysis method [10] to determine the particle charge Q = —(4360 ± 440) e 
and the screening length Xjj = (0.65 ± 0.15) mm. 

At first, without additional heating, the suspension has the solid structure of a triangular 
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lattice with six-fold symmetry. Due to its extreme softness and the stresses applied by 
confining electric fields, this solid is never defect-free. Even ^t the lowest temperatures, it 
has some defects, arranged in strings defining domain walls 



Our laser heating method 



3l| allows adjusting the kinetic temperature of particles by 



varying the laser power. This does not affect the plasma environment or particle charge, 
unlike previous methods 

mm 

. Random kicks are applied through radiation pressure 
from a pair of 532-nm laser beams that are rastered across the suspension in a Lissajous 
pattern with frequences ^ 48.541 H. and /. ^ 30 H. .n a reetangular region sUghtW 

larger than the FOV. During laser heating, the suspension is a driven-dissipative system . 
In steady state, the particle kinetic temperature is determined by a balance of external 
laser heating and frictional drag cooling from neutral gas. Due to the orientation of the 
laser beams, the temperature is higher in the x direction [31] by a ratio of 2 during steady 
heating, and increasing monotonically from 1 to 2 during rapid heating. 

To provide conditions favorable for solid superheating or liquid supercooling, we switch 
the laser on and off abruptly, so that the temperature will change suddenly. In our rapid 
heating and cooling experiment, the pair of 532-nm laser beams is ramped between and 
7 W in 1 or 2 sec, for rapid cooling and heating, respectively. This results in a temperature 
that changes at a rate > 20 000 K/s during rapid heating. 

We measure three indicators of microscopic structure in addition to the temperature time 
series. First, we identify defects and calculate defect area fraction by calculating Voronoi 

pn 

diagrams [3, 130[. Second, we measure short-range translational order using the height of 
the first peak of the pair correlation function g{r) [sil, which is larger for solids than for 
liquids. Third, we measure the short-range orientational order using the bond-angular-order 



parameter, Ge 



321 ], which varies from zero for a gas to unity for a perfect crystal. For a 



solid, Ge is less than unity if there are defects. 

In addition to our measurements with rapid heating and cooling, we also performed 
slower heating and cooling to measure the melting point, in the range 4600 — 5600 K. This 
is consistent with the prediction 4600 ± 1000 K of 2D Yukawa simulations 33|] using our 
measured values of inter-particle spacing, Q and Xd] the error bar arises from uncertainties 
in Q and A^- 

Our main results are time series of temperature and the microscopic structure indicators. 
We applied rapid heating, followed by 55 s of steady conditions and then rapid cooling. The 
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temperature time series, shown in Fig. 1(a), is marked at six times corresponding to the 
Voronoi diagrams in Fig. l(b)-l(g). Time series are presented in Fig. 2 for the structure 
indicators: defect fraction, g{r) peak value, and Gg. We combine time series data to yield a 
hysteresis diagram. Fig. 3. Details are presented next. 

The sequence of Voronoi diagrams. Fig. l(b)-l(g), reveals solid superheating. Before heat- 
ing Fig. 1(b), the suspension has a solid polycrystalline structure, with domains as large as 
the sample region shown here. In the most significant panel in this sequence. Fig. 1(c), at 
T > OOOO-ft' near the end of rapid heating, the structure remains a polycrystalline solid, 
with only a modest increase in defects mostly near the previous defect locations. Since this 
is a solid structure, while at the same time the temperature is above the melting point, we 
conclude that it is a superheated solid. Later, in steady heating. Fig. 1(d), the structure is 
liquid, as indicated by the numerous defects and lack of large crystalline domains. Imme- 
diately after rapid cooling. Fig. 1(e), defects have diminished greatly. Five and ten seconds 
after rapid cooling. Fig. 1(f) and 1(g), respectively, the suspension is again a polycrystalline 
solid, with crystalline domains separated by string-shaped defect clusters. These crystallites 
grow bigger by merging neighbors together gradually in a slow recrystallization process ^. 

Time series for microscopic structure indicators. Fig. 2, reveal different time scales. In 
order to compare these time scales, we rescaled all four variables in Fig. 2(d) to vary from 
(before rapid heating) to 1 (during steady heating) using a linear function with a slope and 
interecept for each variable. During cooling, structure indictors change at different rates: 
translational order changes fastest, and orientational order slowest, consistent with the data 
of 81]. We also measured the defect fraction, which changed at a rate between the other two. 
Our experimental method also allows measurements during rapid heating, where we observe 
a delay in the response of the structure as the temperature increases. This delay is shortest 
(< 0.04 s) for translational order, and longer (^ 0.2 s) for defect fraction and orientational 
order. 

Hysteresis diagrams, like Fig. 3, are traditional tools for studying phase transi- 



tions 
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34l |. Hysteresis arises because structure does not respond immediately to a change 



of temperature. This can occur either due to a delayed response as in the case of our rapid 
heating, or a gradual response as for our rapid cooling. In previous solid superheating ex- 



periments, the vertical axis was typically from X-ray diffraction 



2l|. Here we use direct 



imaging of particles to yield an indicator of microscopic structure for the vertical axis of a 
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hysteresis diagram, Fig. 3. Our hysteresis diagram allows a useful interpretation: a signature 
of solid superheating or liquid supercooling would be a horizontal row of data points across 
the melting point. Such a horizontal row would indicate a temperature that has changed 
without a corresponding change in structure. 

Our chief conclusion, an observation of solid superheating, is based on two results. First, 
Voronoi diagrams compared before and after rapid heating indicate solid superheating, as 
described above. Second, the hysteresis diagram. Fig. 3, has the signature of solid super- 
heating: a nearly horizontal row of data points, which can be seen near the bottom of the 
graph. 

After the superheated solid is formed, it then melts, as indicated by a proliferation of 
defects. Because the substance being melted is a superheated solid, the melting occurs 
without much further temperature increase, yielding a nearly vertical line of data points in 
Fig. 3. The lifetime of the superheated solid and the duration of the subsequent melting are 
both about 0.25 s. 

We also conclude that our rapid cooling did not produce a supercooled liquid. The rapid- 
cooling portion of Fig. 3, lacks the signature of a supercooled liquid. Instead, the defect 
fraction drops dramatically during the temperature decrease. Additionally, the Voronoi 
diagrams for rapid cooling. Fig. l(e)-l(g), lack a liquid structure. 

Our observation that we did not form a supercooled liquid might be attributable to 
the low dimensionality of the experiment. For 3D systems, many examples of materials, 
including dusty plasma [9|], can be quenched to form supercooled hquids or glasses. For 
2D systems, however, forming a supercooled liquid or glass seems to be difficult 25|]. The 
role of dimensionality in transitions to a glassy or supercooled state remains an important 
question [o], [ssl. A previous 2D experiment addressing this question was performed using 
colloidal suspensions 25|], which have much higher friction than in our suspension. 

In addition to our experiment, we also performed a numerical simulation. We found con- 
ditions that result in the same signature of transient solid superheating as in the experiment, 
as we will report in detail elsewhere. 

In conclusion, firstly we have shown that strongly coupled plasmas can exhibit solid su- 
perheating. This suggests investigating superheating in other solid strongly-coupled plasmas 
that can melt, like laser-cooled ions [3i] and the crust of neutron stars f^. Secondly, we have 
demonstrated an experimental method of studying solid superheating using direct imaging 
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of particles. These two results are apparently the first of their kind. Thirdly, we found a 
lack of liquid supercooling in our 2D system. 
This work was supported by NASA and DOE. 
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FIG. 1: (color online), (a) Time series of particle kinetic temperature T{t), when laser heating was 
switched on and then off. Times marked b-g correspond to panels below, (b)-(g) Voronoi diagrams, 
showing defects in color. Polygons indicate the number of nearest neighbors of a particle: red (5), 
white (6), blue (7), and green (others). 



FIG. 2: (color online), (a) Defect fraction (the area of defective polygons, as a fraction of total 
area, in Voronoi diagrams as in Fig. 1). (b) Height of the first peak of the pair correlation function 



32l | , an indicator of short-range 



g{r). This is an indicator of short-range translational order, (c) Gg 
orientational order. This can vary from zero for a gas to unity for a perfect crystal, (d) Time series 
for the three structure indicators and temperature, rescaled to vary from zero (for a solid before 
heating) to unity (for a liquid during heating). Data are smoothed over 3 frames. 



FIG. 3: Hysteresis diagram made by combining data from Fig. 1 and 2. The time interval between 
data points is 0.018 s. Initially, we had a solid, lower left corner. Then rapid heating was applied, 
causing a temperature increase across the melting point without much change in structure, the lower 
horizontal line of data points. This is a signature of solid superheating. Next, the superheated 
solid melted, as shown by the nearly vertical line of data points on the right. The resulting liquid 
in the upper right corner had a high defect fraction. Later, during rapid cooling, defect fraction 
dropped dramatically as the temperature declined. Finally, the suspension slowly recrystallized. 
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